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Abstract 

An integral representation of solutions of the wave equation as a 
superposition of other solutions of this equation is built. The solutions 
from a wide class can be used as building blocks for the representation. 
Considerations are based on mathematical techniques of continuous 
wavelet analysis. The formulas obtained are justified from the point 
of view of distribution theory. A comparison of the results with those 
by G. Kaiser is carried out. Methods of obtaining physical wavelets 
are discussed. 
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1 Introduction 



The aim of the paper is to find a new exact integral representation of 
solutions of the wave equation. We consider a homogeneous equation with 
constant coefficients in a three-dimensional space. We represent a solution 
as a superposition of some other elementary solutions of the wave equation, 
which can be taken from a wide class. Although we use the term 'elementary', 
the solutions may be rather complicated. We call them ^ elementary\ because 
we represent other solutions as their superposition. 

A well-known exact integral representation for the wave equation of such 
a kind is the Fourier integral, where solutions are decomposed into the su- 
perposition of plane waves. However it is sometimes convenient to have 
elementary solutions well-localized in space. This is useful for studying local 
properties of solutions and, among other things, for studying the propaga- 
tion of singularities. An approximate integral representation of solutions as 
a superposition of localized Gaussian beams is developed in However, 
it is heuristic and inexact. Other mathematical methods for developing an 
exact representation of this kind are necessary. 

In our paper, we use mathematical techniques of wavelet analysis. It has 
been developed in the 80s of the 20th century. The main ideas of wavelet 
analysis take their origin in group representation theory and in the theory 
of coherent states (see [51 E] and references therein). The first papers on 
continuous wavelet analysis theory were motivated by applications to seismic 
wave propagation [3 [8]. These papers stimulated interest in wavelet analysis. 
Nowadays a great many of books and articles on wavelet analysis are available 
(see |n]-[Il] for instance). 

Continuous wavelet analysis provides a reconstruction formula, which al- 
lows us to represent functions as a superposition of a family of "elementary" 
functions obtained from one function in a special way. The techniques of con- 
tinuous wavelet analysis imply the following steps. First, we need to fix the 
Hilbert space of functions under consideration. Second, we choose a certain 
function called a 'mother wavelet' in this Hilbert space. The function must 
satisfy a special admissibility condition. Next we specify a group of transfor- 
mations, which will be applied to the mother wavelet in order to construct 
the family of wavelets. In this way, an overcomplete set of functions is ob- 
tained. The group must have some special properties (see [H [lO]). Next 
a wavelet transform of an arbitrary function from the Hilbert space under 
consideration is obtained as the scalar product of the function and each of 
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the wavelets from the family constructed. The wavelet transform depends 
on group parameters and contains information about local properties of the 
function (see [10]). 

The wavelet transforms of all functions from the considered Hilbert space 
themselves form a subspace of the Hilbert space of functions of group param- 
eters. The map from the space of functions to the space of wavelet transforms 
is an isometry. This allows us to obtain a reconstruction formula. We follow 
this scheme in our argument below. 

A special case of continuous wavelet analysis based on the analytic signal 
transform in a three-dimensional space was first applied to the wave equation 
by G. Kaiser in [H] and developed in [I5]-[I7]. He obtained an integral 
representation formula for solutions of the homogeneous wave equation as a 
superposition of elementary solutions derived from one fixed mother wavelet 
only, named by him the 'physical wavelef. 

In outline, the content of the paper is as follows. First we present in 
Section [2] a brief discussion of the idea of our method and illustrate it with a 
few simple facts. A more detailed development of the method is a subject of 
other sections. The aim of Section [3] is to construct a wavelet-based integral 
representation for solutions of the wave equation. We split the whole space 
of solutions into a direct sum of two subspaces containing solutions with 
positive and negative frequencies. In each of them we obtain a family of 
elementary solutions by applying transformations to a certain solution at a 
fixed moment of time. This solution must satisfy an admissibility condition. 
The transformations are spatial translations, dilations, and rotations, which 
are common for standard wavelet analysis. We also apply a dilation of time. 
We suggest to use the wavelet transform of the solution as coefficients for its 
integral representation. We also show that the coefficients can be expressed in 
terms of the solution itself and its time- derivative at a fixed moment of time, 
and this does not require the decomposition of the solution into positive- 
and negative-frequency parts. A justification of the results from the point of 
view of the theory of distributions is given. A detailed comparison with the 
results of Kaiser has been carried out in [12] , and a brief review of it is given 
here. 

In Section m we discuss the possibility of obtaining new physical wavelets 
by means of some known methods of constructing explicit exact solutions of 
the wave equation. We also consider exponentially localized physical wavelets 
found and generahzed in [HI 1201 HI] , wavelet properties of which has been 
studied in [22] ■ 
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2 Preliminary discussion 



Fourier analysis allows one to construct a solution u{r, t) of the wave equation 



uttir,t) -c^[u^^{r,t) +Uyy{r,t)] = 0, 



const, 



(1) 



as the superposition of plane waves. For example, a solution with positive 
frequencies has the decomposition 



u{r,t) 



{2n) 



—- [ d^ku(k) exp(ik ■ r — ilklct). 
T^r J 



(2) 



Here u{k) is the Fourier transform of the function u{r,0). For the sake of 
simplicity, here we consider the two-dimensional case. The aim of our work is 
to construct another integral representation of a solution of the wave equation 
as the superposition of its localized solutions from a wide class. Instead 
of using Fourier analysis, we invoke mathematical techniques of continuous 
wavelet analysis. First we consider a special case of wavelet analysis to 
compare it with the Fourier transform. It allows us to represent an arbitrary 
functions in terms different from harmonic exponents. In the case of the 
Special choice of the mother wavelet, the representation of the function u{r) 
can be written in the form 



uir) 



1 



d^k j d?h Uik, b) exp 



\k\ 



2a^ 



^ifc-(r— 6) 



(3) 



where U{k, b) is defined by the formula 



U{k,b) = \k\ / d ru{r) exp 



2a2 



^--ifc.(r-b) _ g-(T2/2 



, (4) 



and is a constant that plays the role of (27r)2 in the Fourier formula 
(l2l). If 0" ^ 1, the second term can be neglected, and this formula can be 
interpreted as the Fourier transform of the function u{r) multiplied by a 
cutting exponent in the vicinity of the point b. The larger |fc| the smaller 
the width of the cutting function. 

Now we introduce a solution (f'^'^{r,t) of the wave equation, which is 
satisfies the following initial conditions 



\k\ 



2a^ 



\r — b\ 



gifc-(r-5) _ g-a2/2 



(5) 
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|^''V,i) =0. 

Then taking into account formula ([3]), we can obtain a solution u{r, t) of the 
wave equation as follows: 

u{r, t) = ^jd'kj d'b U{k, b) ^'^'''{r, t), (6) 

where U{k, b) is defined by the formula (jl]) and 

u{r,t)\t=o = u{r). (7) 

However we cannot find an explicit solution ^p^'^{r,t) satisfying ([5]). There 
exists an exact highly localized solution of the wave equation named the 
'Gaussian wave packet', which was found in [2U] and studied in [22] and 
which has the following explicit form: 

¥.(r,t) = ^=^==exp(-m/l--), = x - ct + . . (8) 

\/x^ ct-\e \ V e j X + ct - le 

It was shown that if p ^ 1, this solution is actually an exponentially localized 
wave packet moving along the OX axis with speed c. Our aim is to construct 
an integral representation that expresses solutions in terms of the Gaussian 
wave packet ([8]) and other localized solutions of the wave equation. 

The formulas ([3]) - (jl]) are formulas of continuous wavelet analysis but 
written in the nonstandard notation. The notation common for wavelet 
analysis is as follows. The spatial frequency parameter h is viewed as A; = 
(cos/?,sin/?)/a and a,/3 are taken as parameters. The formula for a wavelet 
family reads 

,...V).1.(m-I1^). m,^(-^« -■"/). (9, 

where a defines a dilation, b defines a translation, and /3 defines a rotation 
of the argument of the function yj. Continuous wavelet analysis allows one 
to represent an arbitrary square-integrable function u as the superposition 
of wavelets y?: 

2-K +00 

n{T) = ^ld(3 1^ I d'bU{a,P,b)^^'^'\r), (10) 

IR2 
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M2 



where the coefficients of the decomposition are defined by the formula 

U{a,(3,b) = j d?ru{r) if^'^'^ir). (11) 

The integral for the coefficient C^p must be convergent. This condition re- 
stricts the family of functions what can be used in these formulas. 

Our idea is to take a solution of the wave equation and to use it as a 
wavelet in the initial moment of time. Further the following representation 
gives a solution of the wave equation: 

2tt +00 

u{r,t) = -^1^^ J d'bUia,(3,b)lf^'^''{r,t), 

"^0 ]R2 

uirM=o = u{r), lp'''^'''ir,t)\t=o = v'^'^'^r), (12) 
and the family [p"-'/^'^ is defined by the formula 

^«'/^'V,^) = VfM^^^,-V (13) 
a \ ^ a a J 

The formulas written above are only an illustration of our method. The 
question is whether a solution of any initial- value problem can be represented 
in a similar form. We want to present these formulas in a common wavelet 
analysis formalism. There is also a problem of a convergence of integrals. 



3 Integral representation for solutions of the 
wave equation 

Generally speaking, we seek an integral representation of solutions of the 
wave equation in the form 

u{r,t) = J df,{u) U{u)^%r,t), (14) 

where is a set of parameters, J d//(z/) denotes integration with respect to 
the measure /x(z/) in the space of parameters, (p^{r, t) is a family of elementary 
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solutions dependent on the parameter u, and the (/{u) are coefficients. In 
the following sections, we define each of these objects and the Hilbert space 
7i of solutions where such a representation is allowed. 

3.1 The space Ti of solutions of the wave equation 

Consider the homogeneous wave equation in with a constant coefficient 
c: 

Du = uu - {u^a: + Uyy + u^^) = 0, r = {x,y,z). (15) 

We fix the space 7i of complex-valued solutions of the wave equation as a 
space of functions u{r,t), which are square integrable with respect to the 
spatial coordinate r when the time t is fixed and have the following Fourier 
transform calculated with respect to the spatial coordinates when the time t 
is fixed: 

M(fc,t) =n+(fc,0)e-'l'=l^* + M_(fc,0)e'l'='"*, M±(fc, 0) e L2(M^). (16) 
If the integrals 

j d^k\u±{k,0)\ < oo, J d^k\k\^\u±{k,0)\ < oo, (17) 

converge, the solution u{r,t) allows one to take the second derivative in the 
classical sense and the function u{r,t) is a classical solution of f|T5l) . If the 
integrals f|T7j) diverge, this means that the solution u{r,t) has discontinuities 
in the variable r, and thus is not a classical solution of the wave equation. In 
that case, we introduce solutions in the sense of distributions. The function 
u{r,t) is a solution of (fTS!) in the sense of distributions [23j if it satisfies the 
equation 

^(n(r,t), P{r)) = c'{u{r,t), A(3{r)) (18) 

for all test functions /3(r) in a certain class. We consider test functions 
(3{r) such that they decay at infinity faster than \r\"^, Vm > 0, and have all 
derivatives. This class of test functions is usually named the Schwartz class 
§(M^) (see |25). The notation {u{r,t), P{r)) stands for the common L2(M'^) 
scalar product. If the function u{r,t) G L2(]R^) satisfies the equation f|T8l) 
and have the second derivative in r and t, it is a classical solution of the wave 
equation f[T^ . We note here that even if the function u{r,t) itself cannot 
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be differentiated in t, it is easy to show that the integral {u{r,t), l3{r)) is a 
smooth function of t. 

The space of solutions 7i is decomposed into a direct sum of two subspaces 
7i± defined as follows: 



= : X M ^ C I M+ G n, u+{k,t) = M+(fc,0)exp(-i|fc|ct)} , 

H- = {u_ -.M.^ xR^C \ u_ e n, u^{k,t) = ?i_(fc,0)exp(i|fc|ct)} . 

The spaces T-C± consist of solutions supported on the positive-frequency and 
negative-frequency light cones, respectively. 

In the space we introduce a common L2(M^) scalar product with 
respect to the spatial coordinates r: 



{u+, v+) = / d'ru+ir,t)v+ir,t) = —— / d'ku4k,0)v+ik,0). (20) 



This scalar product does not depend on the time t. This is the main reason 
why we decompose the whole H into a direct sum of these two subspaces 
H±. If we try to use the common L2(M'^) scalar product directly in the space 
Ti, the exponents exp(i|fc|ct) and exp(— i|fc|ct) do not cancel, and the time 
dependance is not removed. The case of 7i_ is analogous to H+. 

3.2 Wavelet-based integral representation for solutions 
from ?{+ 

In this section, we give a decomposition of solutions from 7i+ in terms of ele- 
mentary solutions of the form (IT^ . First we determine a family of elementary 
solutions of the wave equation f|T5|) . 

We fix a solution ip^{r,t) of the wave equation (HM that belongs to the 
space H+. The only constraint here is the following admissibility condition 
on the Fourier transform (p+{k, 0) of ip+: 



We may call this solution a 'physical wavelet', following the terminology 
introduced by G. Kaiser [T?]. We also assume that the physical wavelet 



n = n+®n 



U = Uj^ + u 



(19) 






(21) 
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belongs to L2(M"^) if the time t is fixed. Then the admissibihty 

condition fl2ip holds if the Fourier transform ip+{k, 0) has a root of any order 
at the origin fc = 0. For the sake of simplicity, we assume that the mother 
wavelet has an axial symmetry with respect to the OX axis. A more general 
case will also be considered below. 

We construct a family of elementary solutions in the following way. We 
apply translations by b G M^, dilations by a > 0, and rotations through 
angles 'di, '62 to the spatial coordinates r. We also introduce the dilation 
of the time t by a. We denote by ip"^{r,t) the family of solutions obtained, 
supplying it with the superscript u: 

H^\{r,t) ^ -^ip^ (u-l^J-^, ^) , z/= (a,6,^i,^92), (22) 



a 



where M^^ is defined as follows: 



M^,^2 = sin^^i cos^^i cos^^a -sin^^a • (23) 




In this notation, the three parameters a, di, ^2 have the meaning of a spatial 
frequency vector qr, where \q\ is proportional to 1/a and the angles define its 
direction in spherical coordinate system. 

Now we define the coefficients U{v) of the decomposition ( |T^ . We sug- 
gest to put these coefficients equal to the wavelet transform of the solution 
M+(r,t) G 71+. Then the coefficients C/+(z/) are defined in terms of the scalar 
product of M+ and ip\: 



(27r)3 



d^k n+(r, 0) exp(ifc ■ b) ^+ {aU^l^^k, O). (24) 



Since the scalar product in H± does not depend on t, the coefficients f/+(z^) 
do not depend on time as well. 

Below we will use a fact common to wavelet analysis techniques. The 
wavelet transform of any square integrable functions possesses the isometry 
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property P, [ID] 



(/(r), g{r)) = —J df,{u) F{u) G{u). (25) 

Here /, g are some square integrable functions, and F, G are their wavelet 
transforms. If functions are solutions of the wave equation, they depend on 
the time t parameter. 

We apply (!25l) for two solutions u+(r,t), v+{r,t) taken at one and the 
same fixed moment. We can write 

1 

•p 



v+) = -^ I d/i(z^) U4u) V+iu), Vm+, v+ e (26) 



27r TT oo 



j d/i(z/) = j ddi j d^2 sint92 d^b. (27) 



IR3 



This property extends to any other moment of time, because neither Uj^{v), 
K|_(i/) nor the scalar product of m+ and f + depend on the time t. The isometry 
property implies the reconstruction formula 

u+{r, t) = -^J di^iu) U+{u) (^^(r, t), (28) 



which holds in the weak sense for any G 7^+. The coefficients U^{v) 
do not depend on r and t and thus the formula fl28|) has the meaning of 
a superposition of elementary solutions Lp^^{r,t). Similar arguments can be 
applied to the space 7i_. The weak sense in the last formula means that 
we can take the L2 inner product of m+ and any solution f+ G 7i+ on the 
left-hand side and the inner product of y?!^ and w+ on the right-hand side 
under the sign of integration and obtain equality. This formula can also be 
interpreted from the point of view of distributions of four variables (see [B|) . 

Below we show that the relation (!28|) is valid not only if we take the scalar 
product with a solution f+(r, t), but also if we use a test function /?(r) from 
S(]R^) instead of This means that we can understand the relation fl28p in 
the sense of distributions: 

(n+(r,t),/3(r)) = ^ j d^^) f/+(z/) (^^(r, t), /?(r)), (29) 
where the function I3{r) is an arbitrary test function defined above. 
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The proof of (12^ follows the line of argument given below. First we note 
that the inner product can be written in the form 



(n+(r,t),/3(r))= / n+(r, t) /3(r) 



= (2^/ d'fc?^+(fc,0)e-^-+('=)*5(fc) = («+(r,0), (3.{r,t)), (30) 

where 

P-{r,t) = j d=^fc^(fc)e''=-^+^'^+('=)*. (31) 

The standard isometry formula fl2S]) gives 

(«+(r,0),/?_(r,t)) = dMz/)f/+(z/)54M), (32) 

where the wavelet transform of the function /3_(r,t) can be written as 

i?_(z/,t) = (/3_(r,t), <^:t(r,0)) 



= (/3(r),^:;(r,t)) = (v.;(r,t),/3(r)). (33) 

Substituting ([33D into ([32]), we obtain (|29D. 

Moreover a relation similar to fl29l) but containing the derivatives in time 
is valid: 



where k is any integer. We take the relation (13011 as a starting point and 
differentiate it under the sign of integral k times. For /5 G S we obtain 

AHu,(r,t),,3(r)) „). ^%^\ . (35) 



The isometry property fl25|) results in 



(^+(^-,0), "^'^^^^^^ ) = ^ / d/i(^)t/+(^)5i'H^^,t), (36) 
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where the wavelet transform B_ can be written in the form 

_ d^(/3_(r,t), v.-(r,0)) _ d^(/3(r), ^.-(r, t)) 

dt^ " dt^ • ^'^^^ 

Substituting (P7|) into we obtain This means that if the function 
(f is a solution of (12^ . then the integral is also a solution in the sense of 
distributions. 



3.3 Simplifications and generalizations of the integral 
formula for solutions 

The number of parameters in the family of solutions can be reduced and 
the formulas (1221) - 0281) can be simplified if the mother wavelet if^ is a 
spherically symmetric one. Then the set of parameters reduces to u = (a, b). 
The admissibility condition (|2T1) takes a simpler form: 

oo 

C;=An j dfc ''^+^^'°^'' <oo, k = \k\. (38) 



The family of wavelets reads 

,jV.O^-i,,,(^.i). (39) 

and the representation has the form 

oo 

47r f da 



d'bU^{a,b)vf{r,t). (40) 







Although this particular case of analysis is the simplest one, it is blind to di- 
rectional properties of wave fields, and this may be inconvenient for practical 
purposes. 

The most general case occurs when the mother wavelet ip possesses no 
symmetry at all. Then we should use the full set of three Euler angles '^2 
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and to determine rotations of the mother wavelet in a three-dimensional 
space. The set of parameters now reads v = {a,b,'t}i,-&2,'&3)- The rotation 
matrix now has the form 

M^,^2^3 (41) 

cos t^i — sin -i?! \ / 1 \ / cos — sin 

sin cos i?! I cos ■^2 — sin ( sin -^s cos t^s 

1 / \ sin^^s 008^2 / \ 1 

and the reconstruction formula reads 
u+{r,t) = 1^ j M,j ^2 sint92 j "^^^ j ^ j d'bU+(u) ^l{r,t). (42) 

''^ m3 

Another generalization of the reconstruction formulas is associated with 
the possibility of using different mother wavelets in calculating wavelet trans- 
forms f/+(i^), V+{i') (see [S], [ID])- Let ijj+{r,t), x+('";^) ^ Upon con- 
structing families of wavelets ip'^, x+ similarly to (!22|) . we calculate wavelet 
transforms in the form 

U+{u)^{u+,r+), V+iu)^{v+,xl)- (43) 

The isometry property is valid in this case as well, but instead of the constant 
in the formula (!26|) we must take the constant 



,.^^f,,,MM!^, (44) 



Then we obtain the reconstruction formula 

u+{r,t) = -^ /d/x(z/)t/+(z/)x+(r,t), (45) 

which holds at least in the weak sense and in the sense of distributions. If 
t = 0, then the formula (H5l) is a well-known wavelet reconstruction formula 
(see [9], [I^). The formulas (140|) and (l42l) can be generalized in a similar 
way. 
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3.4 Initial-value problem for the wave equation 

Since each solution of the wave equation from 1-L± can easily be expressed 
in terms of its initial-value problem, it is useful to obtain formulas for the 
solution of the initial-value problem in terms of localized solutions. These 
formulas will allow one to avoid splitting the solution u into its positive- 
frequency and negative-frequency parts m+ and respectively, in order 
to calculate the coefficients of the decomposition. The splitting requires an 
additional operation - taking the Fourier transform of the solution, which is 
not convenient in some cases. Consider the following initial-value problem: 




du 



0, 



(46) 



w{r) e L2(M^), v{r) G Lsi 
We seek a solution of the form 



^^{r,t) = -^ J d^iiu)UM^U^,t) + -^ J dM^)f/-(^)^-(^,i)- (47) 



In this section, we construct the solution formally. We will justify this solu- 
tion in the sense of distributions in the next subsection. 

We need to find the decomposition coefficient f/±(z^). The expression 0241) 
for U± requires m+ and m_ or their Fourier transforms. We obtain u± from 
the initial-value data. To find an idea for this, we use a well-known Fourier 
formula for a solution of this initial-value problem. It reads 



u{r,t) 



1 



2(27r)= 



w(k) rT7^(^) 

ic fe 



exp(ifc ■ r — i|fc|ct) 



+ 



2(27r)3 
This formula yields 



w{k) 



ic\k\ 



:V(k) 



exp(ifc ■ r + i|fc|ct). 



(48) 



w+(fc,0) = - 



w{k) 



iclfc 



■vik) 



u_{k,0) = - 



w{k) 



ic fc 



-v{k) 



(49) 
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Substituting u± from (HUj) into (^^, we obtain 

1 



d^k w{k)if''^{k, 0) 



(60) 



We introduce solutions ip± ^ 'H± of the wave equation in such a way that 

ip±ir, t) = ±-^^^±ir,t), ^±(fc,0) = -ic|fc|V^±(fc,0). (51) 



Famihes of solutions are built in the standard way: 



then 



u = {a,b,di,^2), (52) 



iflir, t) = ±«^^±(^, t), ^l{k, 0) = -i ca |fc| r±{k, 0). (53) 
The formula ( |50l) in the new notation reads 

1 /" ,T .X : TV a 



U±{y) = ^ / w{r) ^l{r, 0) ^ - / dV v{r) r±{r, 0) 



R3 



il^±(z.)T^V±(^), 



(54) 



where the capital letters W± and V± denote the wavelet transforms of w and 
V with respect to wavelets ({)± and %Ij±, respectively: 



(55) 



Now the reconstruction formula of the form (H7|) for the sum of m+ and u. 
reads: 



M(r,t) 



dyu(z/) 



1 



+ — / d/i(z/) 



iH^_(^) + ^^-M 



y.';(r, t) 
V-{r, t). 



(56) 
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3.5 Justification of results 



In this subsection, we show that the formal solution fl56l) can be understood 
at least in the sense of distribution. This means that we understand the wave 
equation P6l) in the sense fl29l) and we state the initial conditions in the form 



/ / n^ «^ ^\ / ( \ w d(M(r,t),/5(r)) 



dt 



{vir),Pir)). (57) 



t=o 



We use the fact that u{r, t) u{r, 0) in the L2 sense as t — > 0. This allows us 
to write {u{r, 0), /?(r)) instead of {u{r, t), (3{r))\t=Q. However we cannot write 
{du{r, t)/dt, P{r)) in the sense of a scalar product in L2, because the function 
du{r,t)/dt in general does not belong to L2 and belongs to the class of 
distributions. Instead, we use the fact that the scalar product {u{r,t), P{r)) 
can be differentiated with respect to t even if the solution u{r,t) itself is 
not differentiable (see \23\). This allows us to pose correctly the initial value 
problem even for nonclassical solutions. If the functions ip± are solutions of 
the wave equation in the sense of (!29|) . the integral in (156|) is also a solution of 
the wave equation in the sense of distributions, as was shown in Section 13.21 
Now we show that u{r,t) defined by (ISBj) also satisfies the initial conditions 

It is useful to take into consideration the auxiliary solution x±(^) t) ^ 

X±(r, t) =±^(^±(r,t), x±(fc,0) = -ic|fc|^±(fc,0). (58) 

The family of solutions x±y constructed by analogy with fl52l) . satisfies the 
relations 

X^(r, t) = ±a^^V^l{r, t), x±(fc, 0) = -i ca |fc| ^^(fc, 0). (59) 

First we consider the simplest case where ip_{r,t) = (p^{r,—t), ip+{r,t) = 
ip{r,t). Then = C+ = C", W+ = W. = W, V+ = V. = V. We fix the 
time t equal to zero and, upon regrouping the terms, obtain 

{u{r,0),/3ir)) = ^^J d/i(^) [Wiu) + Wiu)] (y,'^(r, 0), /?(r)) 
+ / '^^^''^ [-aV{u) + aV{u)] {^^{r, 0),/?(r)) = {w{r), [3{r)) . (60) 
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Here we used the isometry property (1251) for w{r) and /5(r). Next we consider 
the time derivative at the moment of time t = 0. We note that <ff{r, t)|t=o = 
—ip^{r, —t)\t=o. Upon regrouping the terms, we obtain 



^^{u{r,t),(3{r)) 



t=Q 



2a 



d/i(z/) W{u) 



|(v.^(r,t),/?(r)) 



i=0 



t=o. 



2a 



9t 



(^^(r,t),/3(r)) 



t=OJ 



1 



Here we use the fact that 



d/i(z/) V{u){x''{r,0),(3{r)). 



t=o 



t=o 



(61) 



-(¥.^(r,t),/?(r)) 



t=o 



which follows from the Plancherel equality and the fact that dip^/dt = 
is square integrable in r. The last-mentioned integral gives us the function 
(f (r),/?(r)) in accordance with the formulas (HSl) . The constant C^^ can be 
calculated, using (ISTl) and (l58l) : C^x ~ "C*^- 

Now we consider the general case where ip+{r,0) ^ (y9_(r,0). We under- 
stand all the relations below in the sense of ( |29l) and and their analogs 
in Ti.-. In contrast to the previous case, we cannot regroup terms in such a 
way that the expressions under the sign of integration cancel. We must show 
instead that the integrals cancel. First for the time t = the integrals in fl56l) 
multiplied by /5(r) containing terms W^{v)/2 and W^{v)/2 yield the prod- 
uct {w{r), (3{r)) by a wavelet inverse transform formula (the formula fH51) for 
t = 0). We show that other two integrals containing aV^_i_(z/)(y9^(r, 0), /9(r))/2 
and —aV-{i'){ip'^{r, 0), /?(r))/2 cancel and provide the zero sum. Taking into 
account the formulas (159!) . we obtain 



(27r) 



al3{k) 
ica\k\ 



X^(fc,0) = B, 



(62) 
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where Bj^{v) is the wavelet transform of a certain function I3{r) that has the 
Fourier transform of the form f3{k) /ic\k\: 

Mr) = j^,Jd^k ^exp(-ifc.r). (63) 

The singularity at the point = is integrable because the function (]{r) G 
S, and thus j3{k) is continuous. On applying the isometry relation to the 
integral containing aV-^{i>)/2 in (156!) . we obtain 

^ J dM^) (^^(r, 0),/3(r)) = -^(^;(r,0), ^(r)). (64) 

Here we use the fact that C^^ = by the definitions of ip and x in (I5T1) 
and ( l58l) . The integral containing aV^_(z/)/2 gives us the same term but with 
an opposite sign. Their sum is equal to zero; then the formula ( l56l) for t = 
provides (w(r), /5(r)). 

Now we calculate the time derivative of the expression (156|1 multiplied 
by I3{r) at the moment of time t = 0. Here the sum of terms containing 
aV"+(z/)/2 and —aV-{v)/2 yields (t>(r), /3(r)) for the same reasons as those 
concerning the formula ( 16T1) . For example, in the positive-frequency case we 
have 

^ U^i{v)aV^{v)(^y_,{r,t) ,/3(r: 



C+J ^ ^ \dt 



t=o 

d/i(z/) \/+(z/)(x^(r,0),/5(r)) = (t;(r), /3(r)). (65) 
Next we show that the sum of integrals containing W±{v)/2 is equal to zero: 



C+ 



^ U^{u)W^{u)(^yAr,t) 



2C+ / ' ' \ dt 



2C+ 



By analogy with (!62|) . we have 
(xUr,0),/3(r)) 



a(27r) 



d^fcica|fc|^:;(fc,0)/3(fc) 
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j d^fc^!;(fc,0)ic|fc|/3(fc) = B+{v), (66) 



where Bj^{v) is the wavelet transform of a function /3(r) that has the Fourier 
transform of the form (3{k) ic\k\: 

P(r) = J d^k (3{k) ic\k\ exp(-i k-r). (67) 

The integral containing W-^^u) in (!56|) multiplied by j3 is equal to half the 
scalar product {w{r), (3{r)) . The integral containing W-i^u) and multiplied 
by P gives the same term (w(r), /3(r))/2 but with an opposite sign. Then 
their sum is equal to zero and the time derivative of the expression ( !56l) mul- 
tiplied by /3 taken at the moment t = is equal to the function (v(r), /3(r)). 
Then the formula ( |56l) actually yields a solution of the initial-value problem 
fH6|) . Since each solution of the wave equation from H can be represented 
in terms of its initial-value problem data, the integral representation is valid 
for any solution u ETi. 



3.6 Comparison with the results of Kaiser 

We compare our formulas (!22|) - (l42l) with the results obtained by G. Kaiser in 
[Hj. He defines the coefficients of the decomposition in terms of the analytic- 
signal transform (AST) fl9ip of the solution under consideration. Further 
the AST of the solution is regarded in a spatial frequency domain and is 
interpreted as the scalar product fl92|) of this solution and the wavelet fl9^ . 
Therefore the formula for this wavelet is strictly determined by the expression 
for the AST. The parameter cc G of the AST fPTl) has the meaning of 
translation and the parameter s G M, s 7^ 0, has the meaning of dilation. No 
rotation is used, owing to the spherical symmetry of the wavelet (!7i|) . 
The representation formula (1931) obtained by Kaiser then coincides, up to 
notation and normalization, with our special formula PUI) . It should also be 
noted that Kaiser uses another norm fl9Up and do not decompose the whole 
space H. into 1-L± explicitly. 

The main difference between our and Kaiser's approaches is that we start 
our definitions with the decomposition of the space Ti into 7i+ and 7i_ and 
the introduction of time-independent scalar products in each of them, instead 
of using AST. This allows us to choose a mother wavelet from the wide class 
fl^T]) rather than the from fixed one (IMl) used by G. Kaiser |14j . 
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4 Some examples of physical wavelets 



As is clearly seen from Section 13. ![ we can construct the physical wavelet 
for 1-L± just by choosing an arbitrary square integrable function ^±{k) of 
three variables fc G M'^ having a root of any order at the point fc = 0, 
then by multiplying it by the time-depending exponent exp(^i|fc|ct), and by 
taking the Fourier inverse transform with respect to the spatial frequency 
coordinates k. This provides admissible physical wavelets in 7i±. However, 
in practice we possibly will be unable to take the integral analytically and 
find an exact formula for the wavelet in the position space. There are 
several practical methods that allow one to obtain an exact solution of the 
wave equation directly in the position space, without integrating over the 
whole M^. These methods were originated in papers [26]-[30] (see [31j for a 
review of such methods). The aim of this section is to look at some of these 
methods from the point of view of physical wavelets. We find conditions for 
these solutions to be admissible physical wavelets. 

4.1 Spherically symmetric mother wavelets 

G. Kaiser uses in [13] the following method for constructing his physical 
wavelet. He considers two solutions of the inhomogeneous wave equation, 
namely, the emitted wave e(r,t) and the absorbed wave a(r,t): 

(^-c=A)e(r,i) = ^(rt)iM. e(..«) = ji-;^^£^, (68) 

g-c^A)a(r,«)=^(d)i(r), a(r, *) = jlj^^^^^^, (69) 

where is a function of the time t. These solutions are spherically symmetric 
and have singularities at the origin r = 0. The solution e(r, t) has the mean- 
ing of a wave emitted by a point source at the origin, and the solution a(r, t) 
has the meaning of a solution absorbed by a point source. The difference 
between these two functions 

^{r, t) = a(r, t) - e(r, t) = [0(ct + |r|) - 0(ct - |r |)] (70) 

is a solution of the homogeneous wave equation f|T5l) . Upon subtraction, the 
singularity at the origin r = cancels. The function 0(ct) was called by G. 
Kaiser in a ^proxy wavelef. 
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We find conditions that should be applied to the class of proxy wavelets 
in order to obtain admissible physical wavelets. The Fourier transform of 
the solution (p{r,t) flTOjl can be calculated exactly. It reads 



2i|fc|c2 



exp(i|fc|ct) + 



2i|fc|c2 



exp(— i|fc|ct), 



(71) 



where 0(±|fc|) is the Fourier transform of the function 0(t) taken at the point 
±|fc|. This formula shows that the Fourier transform (p{k,t) splits into two 
terms, depending on exp(— i|fc|ct) and exp{i\k\ct). 

We stress here that if we choose a progressive proxy wavelet (f){t), i.e., 
0(^) = for ^ < 0, the wavelet ip constructed will belong to H- only, and it 
can be marked with subscript — , i.e. (p{r,t) = ip_{r,t). The second wavelet 

G ?^+ in that case can be obtained from ip_ by changing the sign of the 
time variable t. The admissibility condition (|38|) then can be stated in terms 
of the proxy wavelet (p in the way 



+ 00 



c- 



dk 



< oo. 



(72) 



The physical wavelet flMl) introduced by Kaiser in [M] was interpreted 
also as a solution derived from the following proxy wavelet: 



Via) 



7r(l -ity 



0(0 = 29(0 r-'exp(-0 



(73) 



where is the Heaviside step function. In the position space, this wavelet 
has the form 



T(a) 



1 



[l-i{ct+\r\ 



[1 - i{ct - \r\)Y 



(74) 



Since the proxy wavelet 0(t) can be chosen from a wide class of functions, we 
can obtain other spherically symmetric physical wavelets with better prop- 
erties from the point of view of wavelet analysis. To show this, we provide 
here another example of spherically symmetric physical wavelet, which has, 
in contrast to that suggested by G. Kaiser, an exponential decay and an 
infinite number of zero moments. This example of a solution first appeared 
in [22], and in [TS] was first regarded as a mother wavelet. It can also be 
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interpreted in terms of a field of two point sources fl68l) . fl69l) with the proxy 
wavelet 

= exp (^-2Vl - it) . (75) 

The branch of the square root with positive real part is implied here and 
below. The difference between the absorbed and emitted waves can be taken 
as a mother wavelet in 7i_ and has the form 



1 



exp 



-2^/l-i{ct+\r\) 



exp 



-2v/l-i(ct- |r|)] } . 



The Fourier transform (p^{k,t) reads 



ifci 



1 

ifci 



(76) 



(77) 



The coefficient C ( |38!) for the wavelet ^p- can be calculated exactly: 



a: 



An J dk 





k 



^5(4) < 00, 



(78) 



where -^'5(4) is a McDonald's function [33] of order 5. The Fourier transform 
of this wavelet has a root of infinite order at the origin fc = owing to 
the factor exp(— 1/|A;|) and the wavelet itself has an infinite number of zero 
moments. The wavelet ip- has a spherical symmetry and an exponential 
decay away from the circle where |r| = ct. 



4.2 Nonsymmetric mother wavelets 

We discuss here the construction of nonsymmetric solutions of the wave equa- 
tion f|T5l) following papers [291 EOj. The method is based on the summation 
of well-known nonstationary Gaussian beams [Ml [35] multiplied by a weight 
function 0(g): 

+00 

^ir,t) = j dq^{q)ipbc^m{q,r,t), (79) 
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exp [iq9{r,t)] 



\Jx + ct — \e\\Jx + ct — ie2 ' 



g > 0, 



(80) 



Q{r,t) = x-ct + 



r 



+ 



X + ct — iei X + ct — ie2 



where ei and 62 are free positive parameters. If 0(g) = 0, g < 0, the formula 
fl79|) has the meaning of a Fourier inverse transform and can be written in a 
simpler form 



1 



^/x + ct — ieiy^X + ct — 182 



(81) 



where is a Fourier inverse transform of (f). The function (p then can be named 
a 'proxy wavelet' following the terminology of G. Kaiser. The formula (|8T|) is 
a special case of the class of solutions presented by Bateman in [261 1221 and 
further developed by Hillion in [28]. Now we determine the class of proxy 
wavelets that produce admissible physical wavelets. The Fourier transform 
of a solution (p{r,t) defined by ( IHTj) reads 



(pik,t) = 271% 



-ilfclct 



k^Ei + kle2 



(82) 



Substituting this expression into the formula for the coefficient C<^ defined by 



(PTj) . we conclude that 0(g) must have a root of order at least 1 + a, a > 0, 
at the origin g = 0. We also restrict the class of admissible proxy wavelet to 
the class Li(M) n^2(K). 

A special case of solutions of the class ( ISTl) named the Gaussian wave 
packet was found in [19], [20] and studied in [22] : 



(^(r, t) 



^Jx + ct — iei\/x + ct — 182 



exp 



1- 



ie{r,t) 



7 



(83) 



where p and 7 are free positive parameters. This solution can be obtained 
from (IHTi) by using the following proxy wavelet: 



(pit) = exp 




(84) 
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The Fourier transform of the Gaussian packet 
the form 



due to 



and 



has 



X exp 



-1 



p 



^7 |fc|(|fc| + fcJ3/2 
2 1 



2 ' 27|fc| + A;x 2{k^+\k\) ' ^^^"^ 

This physical wavelet has an exponential decay away from the moving point 
X = ct,y = 0, z = 0. It has infinitely many zero moments with respect to 
spatial coordinates. As is shown in [22], its asymptotics coincides with the 
Morlet wavelet [21 ED] as p — >• oo and the time t is fixed: 



C 



exp 



;-i£l)V2(_ie2)l/2 

X [1 + 0{p-^"+^)] 



i/t(x — ct) 



[x 



ctf 



2^^ 

a G (1/3, 1/2), 



where 



47 Vp, (^l = 7^1 /P, (^1 = 7^2 /P, 



2al 



P_ 

27' 



2(t2 



(86) 
(87) 



in the domain 

(a; - ct)/7 = O (p-") , y/^=0{p--), z/^=0{-p^), 

provided that the parameters Ictjej, p~^'~^/ej, j = 1,2, are small. 

The axially symmetric case of the Gaussian beam (I8OII and the Gaussian 
packet (|83|) can be obtained by putting ei = 62- 



5 Conclusions 

A new integral representation of solutions of the wave equation was built. It 
is based on mathematical methods of a continuous wavelet transform in a 
three-dimensional space. An arbitrary solution of the wave equation can be 
represented as a superposition of elementary solutions. We discussed meth- 
ods of constructing these elementary solutions. Both spherically symmetric 
and axially symmetric elementary solutions were considered. Examples of 
known physical wavelets were given. Solution of the initial-value problem 
based on wavelet analysis was also considered. A brief comparison of the 
presented results with the results obtained by G. Kaiser was also carried out. 
The integral representation constructed in this article may be useful in real 
problems of wave propagation. 
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A The representation build by G. Kaiser 

The first integral representation of the form (1141) based on mathematical 
methods of wavelet analysis was presented by G. Kaiser in his book [H]. We 
give here a brief review of his results for the scalar wave equation (or acoustic 
equation). The results for vector Maxwell equations were also presented in 
|14j . In this Appendix we follow the specific notation introduced by Kaiser. 

A solution of the wave equation can be represented in the following form 
by means of the Fourier transform: 



x = {x,t), £c G R^, P = ip,Po), p G M^, uj = \p\. 

Here, in accordance with the notation by G. Kaiser, F{x), x = {x, t) denotes 
a solution of the wave equation and /(p), p = (p, Po); is the Fourier transform 
of F{x) taken with respect to both spatial and time coordinates, using the 
Lorentz-covariant scalar product p ■ x = p^t — p ■ x. Owing to the wave 
equation, the function /(p) depends not on four variables p,po but on three 
variables p in the following two ways: f{p, \p\) or f{p, —\p\)- 
The norm of the solution F is defined in the way 



The analytic signal transform of F{x) plays the role of coefficients for the 
decomposition f[T^ : 




(89) 




(90) 



oo 





— oo 
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Kaiser showed that it is sufficient to put z = x + iy equal to (a;, is), x e 
M^, s 7^ 0. Then the analytic signal transform (1911) takes a simpler form in 
the Fourier domain: 

F(x,is)= / /!^c."-ie-P- [e(s)e--V(p,c^) + 0(-^)e"Vb,-^)] ,(92) 

where G is the Heaviside step function. Further an arbitrary solution -F(x) 
with the finite norm fl90|) can be represented in the form 

F(x') = j d/i«(z) ^.(x') F(z), x' = {x\ t), z = {x, is), (93) 



E 



d/ia(z) = T^T^-A^Td'a^dslsr-^, E = {(a;,is), a; G M^ s ^ 0}, 
I \a — Z) 

^a3is(x')= /" ^i^u;"-ie'P-(^-^') re(cus)e-'^^+'"* + 0(-cjs)e'^^-''^*l . (94) 
J Stx-^uj 

The integral (1931) is a special case of the representation in the most general 
form (I14p . The coefficients -F(z) do not depend on time and coordinates, and 
thus the formula (p3|) has the meaning of a superposition of solutions 
The solution is also considered in Section 14.11 in connection with proxy 
wavelets. In the position space it is given by the formula (|74l) . 

A disadvantage of the formula (|93l) is that we can use only one solution 
\E'(x') as a 'building block' for the reconstruction. The method presented in 
this paper is free of this disadvantage. 



B Relationship between the scalar product in 
and four-dimensional distributions 

In Section 13.11 we introduced a solution of the wave equation in the sense 
of distributions. We regarded solutions as distributions in three spatial vari- 
ables r = (x, y, z) and the time t that is viewed as the parameter. This 
approach is similar to that presented in [23] • However we may view solu- 
tions as distributions in all variables r,t. Then we must use four-variable 
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test functions a{r,t) G S(]R'*). The Fourier transform of a solution u & Ti 
calculated in four variables has the form 

F[u]{k,uj) = 5{uj + c\k\) u+{k,0) + 5{uj-c\k\) u.{k,0). (95) 

Next we calculate the action of u on a{r,t) G S(]R^), using the Plancherel 
equality {u{r,t), a{r,t)) = l/27i^ {u{k,Lj), a{k,Lj)). We obtain 



u{r,t),a{r,t)) = ( M+(fc,0), a{k,c\k\) ) + ( M_(fc,0), a{k,-c\k\) ) . (96) 



We introduce solutions A^{r, t) and A_(r, t) by means of the inverse Fourier 
transforms of a{k,c\k\) and a{k, —c\k\), respectively, in the form 



^+ = d'ka{k,c\k\)e"'-^-'^\ A. = ^ j d'k a{k, -c\k\) e"'-^^'-'. 

Finally, formula (196!) can be written in terms of these solutions: 

{u{r, t), a{r, t)) = {u+{r, t), A+(r, t)) + (n_(r, t), A^r, t)). (97) 

We emphasize that a{r, t) is not a solution of the wave equation, but A+{r, t), 
A_{r, t), which are constructed by means of a, are solutions. So we conclude 
that the action of u on a G S(R^) can be expressed in terms of the scalar 
product in H. 
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